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We study how the addition of defects in an otherwise uniform Josehpson junction array modifies 
the critical current of the array, by numerically solving the system of coupled Josephson equations. 
Our results confirm the existence of two different regimes, depending on the normalized maximum 
Josephson current of each junction, i max : for large z max , the defects deteriorate partially the critical 
state of the uniform array and so decrease the critical current, while for low values of i max , the 
defects help the formation of the critical state, which leads to an enhancement of critical current. 
The results are discussed in comparison with the behavior of type-II superconductors, for which the 
inclusion of defects always results in an enhancement of the critical current. 



The traditional way of improving the critical current 
of hard type-II superconductors has been the introduc- 
tion of defects, which act as pinning sites of the vortices. 
Based on the analogy between the magnetic properties of 
superconductors and Josephson junction (JJ) arrays, it 
seems reasonable to expect that in the latter system the 
presence of defects will yield to an improved total critical 
current. This has been the predominant point of view, 
starting from the work of Josephson himself . His point 
of view was followed in later researches on defects con- 
taining JJ's ^-||, and in systematic works on planar J J 
arrays, where Josephson vortices (JV's) and their depin- 
ning were studied (for example, see f8-[l2|). Moreover, it 
has become evident that one of the limiting factors for 
the critical current of actual high-T c superconductors is 
the presence of a granular structure consisting of super- 
conducting grains connected by JJ's. The role played by 
the defects in JJ arrays with respect to transport current 
is therefore a topic of practical importance as well. Clem 
introduced the JV pinning mechanism as the reason for 
the intergranular critical state (CS) |13|| . In this work, 
we study this problem by numerically solving the equa- 
tions describing the behavior of a JJ array in which de- 
fects are introduced. Our results will demonstrate that, 
in contrast with the case of type-II superconductors, the 
total critical current in the array decreases when adding 
defects, except in the regime of low Josephson current, 
for which the inclusion of defects may help the formation 
of a CS which does not originally appear in the uniform 
JJ array. 

In order to avoid complicated demagnetizing effects, 
the studied JJ array consists of a number of infinitely 
long superconducting grains along the z axis. Its xy- 
plane cross-section forms an N x x N y square lattice of 
parameter a, each grain being centered at with i = 

1,2, . . . , N x and j — 1,2,..., N y . Every pair of nearest 
grains is weakly linked by a (when looking at it along the 
z axis) short J J, which has a line density of the maximum 
DC Josephson current / max (A/m) and a line specific 



resistance R (fhn), the line being along the z dimension. 
The capacitance effects are neglected for studying a slow 
process. Four nearest grains with an enclosed void form 
a square cell, the effective void area being A v (m 2 ). We 
name the gauge-invariant phase differences (GIPD's) for 
JJ's along the x and y axes 9fj and 6\p respectively. 

When a transport current with line density /array is 
uniformly applied in the x direction, or a field H is uni- 
formly applied in the z direction, a set of differential 
equations for the 2N x N y - N x - N y GIPD's of all JJ's 
of the JJ array is built up based on the DC and AC 
Josephson equations and the Ampere and the Ohm laws 
as 
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In these equations, t* is the normalized time t to the nom- 
inal time constant r of one cell, t* = t/r = tR/fioA v , and 
hj and i max are the boundary field Hj and I max normal- 
ized to $o/mo^v- Correspondingly, the normalization of 
-^array into i ar ray can be made in the same way as I max . 
We have hj = (j/N y — l/2)i aiia , y in the transport case and 
hj = h in the magnetic case. The field Hij in the 
void (i and j are defined as those for the left-and-down 
adjacent grain) produced by DC Josephson currents nor- 
malized to $0/ H$A V is calculated by 

% = - + n,j+i - 0ij)/ 2 ^ ( 2 ) 

where we have assumed the phase of the order parameter 
to be zero in each grain. 

It is known |l4|,[l5| that a uniform slablike JJ array 
jl6| can be in the CS when i max is large enough, while for 
lower values of t max it presents a JV state. Such an evolu- 
tion in ID JJ array from the JV state to the CS has been 
studied mathematically by Greene in terms of a standard 
mapping technique. His conclusion is that a stochastic 
transition occurs at i max = i max = 0.9716 • • ■ /2ir ms 0.155 

The transition from a JV state to a CS occurs in our 
(mathematically 2D) JJ array as well. To show this 
transition, we analyze the field profiles in a JJ array of 
N x = N y = 25 in the saturated remanence state. The 
profiles are calculated using a Runge-Kutta method as 
follows. We start from {9*j v — 0}, step h from to 100 
and calculate the fully relaxed {hij}, and then step h 
down to and calculate again the fully relaxed {h^}, 
which is our solution. 

As seen in Fig. 1, the profiles of uniform J J arrays with 
*max = 0.01, 0.05, and 0.1 show JV states with 8, 24, and 
64 JV's of roughly equal heights. For i max — 0.17, there 
is a central peak remarkably higher than the others, in- 
dicating a transition from the JV state to the CS. When 
imax > 0.4, the entire profile appears to be a single peak 
of pyramid shape (as shown in the figure for i max = 1), 
which is the distinctive characteristic of the Bean CS, as 
in hard superconductors, where the penetrated supercur- 
rents flow with a density J equal to the critical-current 
density J c @. 

Let us now study the effect of introducing defects in 
the array by changing i max of JJ's connecting certain de- 
fect grains into i maX! d = ci max , where < c ^ 1. First 
we consider the case for low values of i max , for which the 
uniform JJ array does not show a CS. In order to be able 
to check the validity of the solutions by their symme- 
try, we choose symmetrically and uniformly distributed 
defects shown in Fig. 2(a) with c = 0.1. The trend ob- 
served in the calculated results shown in Fig. 1 is clear: 
defects help the formation of the CS. For the nonuniform 
JJ array of i m ax = 0.01, there is a 12-peak structure with 
eight outer peaks lower than the inner four. When i max 
is increased to 0.05, the profile appears to be a single 
round peak, indicating a transition from the JV state to 
the CS. The profiles for i max = 0.1 and 0.17 are already 



of an overall pyramid shape, typical for the Bean CS, in 
clear contrast to the results for the uniform JJ array. We 
confirm therefore that the inclusion of defects in the ar- 
ray with low values of i max results in the enhancement 
and even induction of the CS. 

After this conclusion, an important question arises as: 
Is this behavior general for all values of i max ? The pre- 
vious result of the decreasing of the trapped flux for 
«max = 1 [as seen in Fig. l(i)] already gives us a hint 
of the actual striking behavior: for such values of i max 
that the uniform array develops a good CS, the addition 
of defects results in a decrease of critical current. 

To further confirm the above statement we set two 
square defect clusters shown in Fig. 2(b). Our aim is 
to calculate the line density of the critical current of the 
J J array, i c , a rray, namely, the maximum i array at which a 
static solution can be obtained. The reference is chosen 
as the uniform (c = 1) J J array of i max = 1. We start 
from a zero- field cooled state at t* = with {6*- v = 0}, 
and change i ari ay between and 25 (the maximum pos- 
sible value for i max = 1) to calculate the fully relaxed 
GIPD's that corresponds to i c ,array 

The calculated i Cj array as a function of c is given in 
Fig. 3. We can see that the presence of defects makes 
ic, array smaller in general. It equals 24 when 0.995 < c < 
1.519, and decreases stepwise if c is out of this range. 
The reasons for this are twofold, (i) Since i array flows in 
the horizontal x direction through JJ's of i = 1 to 25 in 
series, it is limited by the ideal critical current J c , array, ideal 
defined as the smallest sum of « m ax and i m ax,d of all the 
horizontal JJ's for each ith vertical cross-section. For 
the array in Fig. 2(b), i c ,array,ideal = 25 if c > 1 and 
i c , array, ideal = 19 + 6c if c < 1. (ii) The phase coupling 

related to Eq. (2) makes i c ,array < ic, array, ideal, so that 

the calculated i c . array < 25 if c > 1 and i c , array < 19 + 6c 
if c < 1. The effect of the defects on the critical current 
of the nonuniform array can be observed in greater detail 
from the calculated current distributions in the array for 
different values of c presented in Q . 

To sum up, we have shown that the addition of defects 
in a JJ array helps to settle a CS in the system if the 
normalized JJ current i m ax is small. However if i max is 
large enough, that is for the cases in which the uniform 
array is in a good CS, the introduction of defects deteri- 
orate the critical current of the array. The latter feature 
is in contrast with the popular view that the CS arises 
from JV pinning. From these results, it is imperative to 
conclude that the origin of the CS in JJ arrays is a phe- 
nomenon more complex than a JV pinning analogous to 
the Abrikosov vortex pinning in hard type-II supercon- 
ductors. 

We thank Ministerio de Ciencia y Tecnologfa project 
number BFM2000-0001, and CIRIT project number 
1999SGR00340 for financial support. 



2 



[1] B. D. Josephson, Rev. Mod. Phys. 36, 216 (1964). 

[2] T. Yamashita, L. Rinderer, K. Nakajima, and Y. On- 

odera, J. Low Temp. Phys. 17, 191 (1974) . 
[3] V. M. Vinokur and A. E. Koshelev, Sov. Phys. JETP 70, 

547 (1990). 

[4] R. Fehrenbacher, V. B. Geshkenbein, and G. Blatter, 

Phys. Rev. B 45, 5450 (1992). 
[5] M. A. Itzler and M. Tinkham, Phys. Rev. B 51, 435 

(1995). 

[6] S. Teitel and C. Jayaprakash, Phys. Rev. B 27, 598 

(1983); Phys. Rev. Lett. 51, 1999 (1983). 
[7] C. J. Lobb, D. W. Abraham, and M. Tinkham, Phys. 

Rev. B 27, 150 (1983). 
[8] D. Stroud and S. Kivelson, Phys. Rev. B 35, 3478 (1987). 
[9] W. Xia and P. L. Leath, Phys. Rev. Lett. 63, 1428 (1989). 
[10] M. B. Cohn, M. S. Rzchowski, S. P. Benz, and C. J. Lobb, 

Phys. Rev. B 43, 12823 (1991). 
[11] Y.-H. Li and S. Teitel, Phys. Rev. Lett. 67, 2894 (1991). 
[12] E. K. F. Dang, G. J. Halasz, and B. L. Gyorffy, Phys. 

Rev. B 46, 8353 (1992). 
[13] J. R. Clem, Physica C 153-155, 50 (1988). 
[14] J. M. Greene, J. Math. Phys. 20, 1183 (1979). 
[15] S. N. Dorogovtsev and A. N. Samukhin, Europhys. Lett. 
25, 693 (1994). 

[16] The slablike is defined as a JJ array studied in the present 
work with N v = oo, so each row of grains or JJ's should 
be in the same state when the JJ array is placed in a 
uniform field in the z direction; it is the analogous to a 
superconducting slab. 

[17] C. P. Bean, Phys. Rev. Lett. 8, 250 (1962). 

[18] D.-X. Chen, J. J. Moreno, A. Hernando, and A. Sanchez, 
in Studies in High Temperature Superconductors, edited 
by V. Narlikar and F. M. Araujo-Moreira, Vol. 40, p. 1 
(Nova Science, New York, 2002). 

FIG. 1. Trapped field profiles of 25 x 25 square-columnar 
J J arrays at saturated remanence. (a), (c), (e), (g), and (i) 
are for uniform JJ arrays; (b), (d), (f), (h), and (j) are for 
nonuniform J J arrays with i max ,d = «max/10. Relative h scale 
factors 1/imax are used. 

FIG. 2. The xy cross-section of the studied JJ arrays. The 
defect grains are marked by shadow, (a) and (b) are chosen 
for the cases of small and large values of i ma x, respectively. 



FIG. 3. The JJ critical current i c , array as a function of c. 
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